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Subject 


Abstract 


An  upper  boun^yis  established  for  the  magnitude  of  the  truncation 


error  Incurred  wharoi  a read^valued,  finite  energy  signal  which  Is  band 


itr)(o  <r<l)  is  approximated  by  2N+1  terms  from  Its 


limited  to 


Shannon  sampling  series  expansion.  The  sampling  expansion  is  associated 
with  the  ban<^  J]-it^n^\and  consequently  involves  samples  taken  at  the 


The  bound  obtained  Is  of  the  form 


where  E Is  the  signal  energy.  This  bound  is  of  the  same  asymptotic 


form  as  the  bounds  derived  by  Yao  and  Thomas  and  Brown 


bound  derived  here  is  tighter  than  the  Yao-Thomas  bound  for  values  of  r 


near  unity,  and  Is  tighter  than  the  bound  obtained  by  Brown  for  all  value 


The  attached  paper  Is  being  submitted  for  publication  In  the  IEEE 


Transactions  on  Information  Theory 


I.  Introduction 

It  is  well  known  that  a band-limited  function  having  finite 
energy  can  be  represented  exactly  for  all  times  by  an  Infinite 
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series  involving  samples  of  the  function.  Several  authors 


CO. 


|J2J  have  obtained  bounds  on  the  truncation  error  that  occurs  when 
the  function  Is  approximated  by  a finite  number  of  terms  from 
the  sampling  theorem  expansion,  rather  than  the  infinite  series. 

In  this  paper  another  bound  is  derived  which  represents  an  improve- 
ment over  the  previous  bounds. 

In  the  analysis  that  follows,  we  will  consider  real-valued,  band- 
limited  functions  with  finite  energy.  Without  loss  of  generality, 
we  can  assume  that  the  function  is  band-limited  to  ( -*,*).  That 
is,  the  function  has  a representation  of  the  form. 


p*  lent 

*2*  \ *(©)«  doo  (-00  <t<  oo) 


where  F(o>)  is  a comp lex- valued  function  having  the  property  that 


f I I2  f 

|F(®)|  dcD  - 2*  J 


f (t)dt  - 2*E  < oo. 


By  the  Shannon  sampling  theorem,  f(t)  also  has  the  representation 


oo 

■I 


sinx(t-n) 


*(t-n) 


(-oo<t<oo) . 
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The  problem  that  will  concern  us  here  Is  to  estimate  the 


truncation  error  Incurred  when  f(t)  is  represented  by  a finite 


number  of  terms,  rather  than  the  Infinite  series  given  by  (3) 


We  will  define  the  truncation  error  by 


3 I that  for  this  class  of  functions,  the  magnitude 


It  is  known 


of  the  truncation  error  is  bounded  by 


Thus,  the  truncation  error  goes  to  zero  asymptotically  like  N 


examined  the  bounds  on  the  truncation  error 


Yao  and  Thomas  I 1 


when  the  additional  assumption  was  made  that  F(u>)  vanishes  in 


(-«,  -xr)  and  (xr,Ji),  where  o < r < 1.  That  is,  they  considered 


the  case  where 


and  the  truncation  error  is  defined  by  (4) , involving  samples 


taken  at  the  Integer  points.  Hence,  the  signal  is  sampled  faster 


than  required  for  reconstruction  by  the  sampling  theorem.  Using  a 


contour  Integral  to  represent  the  truncation  error  and  then  applying 


known  results  concerning  the  growth  of  the  entire  function  f(z)  in 


the  complex  z-plane,  Yao  and  Thomas  obtained  a bound  of  the  form 


-4- 
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4 <*> 


Thus,  by  Introducing  the  guardband  (or  equivalently,  by  sampling 
in  excess  of  the  Shannon  rate) , Yao  and  Thomas  obtained  a bound 
which  behaves  asymptotically  like  N * for  large  N. 

Brown  |V|  also  examined  this  class  of  functions,  and  using  only 
results  from  real  variable  theory  , obtained  a bound  with  the  same 
asymptotic  behavior  as  the  Yao-Thomas  bound.  The  bound  obtained 
by  Brown  is  an  improvement  on  the  Yao-Thomas  bound  for  values  of  r 

near  unity;  that  is,  for  samples  taken  nearly  at  the  Shannon  rate. 

1 

In  this  report,  a bound  is  obtained  which  is  an  improvement  on 
the  bound  obtained  by  Brown  for  all  values  of  r,  the  guard-band  coefficient. 

1 

Like  the  analysis  used  by  Brown,  only  real  variable  results  are  used, 
and  no  application  is  made  of  the  rather  deep  results  concerning  the 
theory  of  entire  functions. 


II  Analysis 

Consider  a real-valued  function,  with  a representation  given 
by  (1).  Defining  the  truncation  error  by  (4)  and  using  the  represen- 
tation for  f(t)  given  in  (3),  we  have 


- X *«  ^f1 

-00 


slnw(t-n) 

Jt(t-n) 


X”-  slnn(t-n) 

2_  £<n)  »(t-n) 

-N 


(6) 


Note  that  for  t-k,  where  k is  an  integer. 
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slmt  (k-n) 
n(k-n) 


In  what  follows,  we  will  assume 


< N,  and  that  t is  not  equal  to 


an  integer.  From  (6)  we  have  that 


simtt 


simtt 


Using  the  Schwarz  inequality 


sin  it  t 


sinrtt 


I 


-6- 
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11  ■ ■ 1 1 


W jnjiipmjiiH". 


Substituting  (8)  and  (9)  into  (7) , we  have 


< Isinntl 


[2NE 


1/2 


(10) 


This  result,  which  is  a slight  generalization  of  a result 

derived  earlier  by  Balahrishman  &]  , provides  a bound  for  the 

magnitude  of  the  truncation  error,  in  terms  of  the  energy  of  the 

band-limited  signal,  which  goes  to  zero  asymptotically  like 
-1/2 

N as  N goes  to  infinity. 

We  now  make  the  additional  assumption  that  f(t)  is  band- 
limited  to  (-nr,  nr)  with  o<r<l.  Then,  f(t)  has  the  representa- 
tion given  by  (5) . Defining  the  truncation  error  by  (4)  and  using 
the  representation  for  f(t)  given  by  (5),  we  find 


-7-  , 
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But 


1_ 

2n 


-x 


sln(N+2)  (u-  at ) i cd  t d e 

.m(^)  2" 


1 CD  t 


■ e 


icDt  -ic 


sln(N»2)  (u-  g>  ) du 


(13) 


Substituting  (12)  and  (13)  into  (11)  we  obtain 

?r  * 1 r 


•fc>  ■ h 


0 


-xr 


F<a>)  I i-  f SlA(^.(u^CP)  lent  iut 

v ^ •!»(—)  L 


du  j>  dw  (14) 

J 


We  will  now  restrict  our  attention  to  t in  the  range 
H_-N,  nJ  . As  noted  earlier,  the  truncation  error  is  zero  for  t 
equal  to  an  integer  and  thus,  without  loss  of  generality,  we  can 
assume  that 

t - m + b (15) 

where  m is  an  Integer  such  that 

- N < m < N (16) 

and 

o < |b|  < \ 

Then, 


2 . 


(17) 


i to  t iut  i(nrt-b)cD  i(mfb)u 
e -e  ««  e - e 

i(nri-b)  cd 


ei(mfb)u 

£l  . eim(u-cD)^^(p3ij2j  . 


But 


[H*pO]2- 


r b(u-cD) 

\ cos 2 + 1 8in 


b(u 


H 


- 1-2  .in2  21  .in 


-8-  . 
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Thus, 


i<Dt  -eiut  +2elm^a“m^sin2  b(u*tp) 

-i  8ln  bfcsJL  cos 


Substituting  (18)  into  (14)  and  making  the  change  of  variable  y ■ — — 
in  the  integration  over  u,  ve  obtain 

e^(t)  - ^ | F(  ® ) + I2  - i I3J  du>  , (19) 


where 


it-  Qp 

If2 


-CEhflL) 

2 


sin  y 


(l  - «12"y)  dy. 


sin  y 


-(*H-®  ) 
2 


12my  ,2 
e sin  by  dy 


and  I„ 


sin  y 


-(«+  ® . 
2 


ei2my  J)y  COS  Jjy  <ly# 


Note  that 


X1  “ 2 


sin  y 


( 1 - e12”y)  d. 


i P» 


sin  y 


l-e120*  dy. 


-9- 
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I 


k,  ^ 

§ ! 


■ 


since  m is  an  integer  and  the  integrand  is  periodic  with  period  «.  Also, 
sin  (2N+l)y 


sin  y 


it/2 


is  an  even  function  and  hence, 

f 2 .in  ( . cos  2m  ) A 

1 J sin  y v 

o 

n it 

11  /~2  sin  [2(N+m)+ll  y . f2  sin  [2(N-m)+l]  y , 

2 ( J sin  y * J sin  y * 


sin  (2N+l)y  ^ 
sin  y 


0. 


The  problem  that  concerns  us  now  is  to  evaluate  Ig  and  I^.  The 

approach  used  will  be  to  show  that  for  o<  | b | < -nr  <c»<nr  and  o<r<l, 

sin^  y . . . f ,n+o\  n-uTI  sin  by  cos  by 

is  monotonic  on  I g-)  , — and  - 


sin 


sm  y 


is  monotonic 


r ,n+oi.  ! 

, , r n-05 

[-hr><  °J 

and  also  on  0,  — g- 

integrals  is  then  used  to  evaluate  Ig  and  I^. 


w-  00 
2 


f ‘ ei2m^  sin2  by  dy 

I sin  y J 

J rr 


2 } 


(23) 


-10- , 
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We  are  Interested  in  u>  in  the  interval  £-nr,  Krl . Hence, 


«+t» 


-*  < (1+r)  < - (^-)  < (1-r)  <0 


and 


0 < | (1-r)  < SZSL  < § (1+r)  < «. 

sin^  bv 

But  it  is  shown  in  Appendix  I that  — ■ ■ is  monotone  increasing  on 

2 sin  y 

(-«,«).  Hence,  S-^n  is  monotone  increasing  on  jy-(-^-°),  "J  > and 

by  the  second  mean  value  theorem  for  integrals  (4) , 


I 


2 


. 2.  ,n+(D  s 
sin  b(— ^ — ) 

. ,Jt+'Dx 
sin  (~2 — ) 


sin  (2N+l)y  cos  2my  dy 


o - HUP 

. 2,  #ir-  as  x _—t- 

sin  b(— ^ — ) r 2 

+ j sin(2N+l)y  cos  2my  dy 

sin(^)  J 

Pi 


-i 


, 2,  ,n+<p  s 
sin  b(— 2 — ) 

sin(— j — ) 


sin(2N+l)y  sin  2my  dy 


, 2 £-cd  r 2 
sin  b(— — ) J 

i /Hr®. \ J 

.ln(— ) J 


* — 


1974 

R-1018 


+i 


sin(2N+l)y  sin  2my  dy 


(24) 


where 


Also, 


- (^p-)  < Pj,  < 


,ft+OK  n-  a 

( 2 > - P1  - 2 ‘ 


2 ' 


sin(2N+l) 
sin  y 


ei2my  g£n  by  cos  dy 


f 

J(^> 


sin(2N+l)y  cos  2my  (~ 1^^)  dy 


sin(2N+l)y  cos  2my  dy 


J ,fl+cn . 
2 1 


sin(2N+l)y  sin  2my  ("2sin"y)  dy 
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sin(2N+l)  y sin  2my  ("gVin^'y*  dy  * 


r 


-13- 


April  18, 
HSP: pgs 
Serial  No. 


and 


- (^r6-)  < p5  < » . 


The  integrals  appearing  in  (24)  and  (25)  are  evaluated  in  Appendix  II, 
and  we  obtain 


|T  <T  I (1+  \l~2 ) I sin  bn|  

lI2  " iI3l  - „ 03  l 2N+2m+l  ' 2N-2m+l 


1 + - 1 


2 cos 


(Isin  btn  I + cos  bto  ) 


Eo 


From  (19)  and  (23) , 

eN  (t)  “ ~2  I F(cD)  ei(mfb)(B  (I2  -i  I3)  do>  . 

Jt 

•itr 

Applying  the  Schwarz  inequality, 


■ nr 

nr 

le*  (t)l  < ~2 

[/ 

If(co)  P dco 

f |l2  - i I3I  2 dto 

-i  1/2 


■Jtr 


■nr 


But 


r «r  | 1 

/ \ 

1 ll2  - i I3  1 2 dco  < 

(1+/5)  sin  bn  [ 1 1 

2 \2N+2m+l  T 2N-2mfll 

J 

-nr 

_ \ /_ 

nr 

1 


(Isin  b oo  1 + Icos  b co  I ) ^ 

2 to 

cos  — 


to 


Using  the  Schwarz  inequality  for  sums. 


2 I 2 I j 2 

( Isin  bool  + Icos  bco  I ) < ( sin  bail  + Icos  bio  I ) (1+1) 


197U 

R-1018 


(26) 


(27) 


(28) 


2 , for  all  real  to 


Hence, 


and  thus 


-14-  , 


I I 2 

|sln  bg)  4-  cosbo3  ) do3  < 2 

2 ^ 
cos  _o3_ 

2 


2 03  . 

sec  — d 03 


» nr 
- 8 tan  — j , 


j |l2  - il3|  do3  < 2 ^(l+<fF 


April  18,  1974 
HSPrpgs 

Serial  No.  R-1018 


sin  bn  | 2N+2nH-l  + 2N 


3xi)[c,n  -2  )l/2] 


Also, 


I F(  03  ) I do3  ■ 2nE. 


Substituting  (29)  and  (30)  into  (27),  we  obtain 

1/2  / 


4 <4 


nr 

fnEtan  2 


[1+J2)  sin  bn 


N+nrt-1/2  N-mfl/2 


Recall  that  t •»  mfb,  and  note  that  sinbn  » fein(ra-t) n I = Isin  nt  j, 


14, 


so  that 


k (t>l 


(nEtan  2 ) (1+\|2)  Isin  nt 


N+m+1/2  N-nH-1  / 2 


where  m is  the  nearest  Integer  to  t , and  It  I < N. 

Equation  (31)  is  our  desired  result.  This  bound,  involving  the 
square-root  of  the  signal  energy,  has  the  same  asymptotic  behavior  as 
the  bounds  derived  earlier  by  Yao  and  Thomas  M and  Brown  [2]  . That 
is,  all  three  bounds  are  of  the  form 


l 


" ■'  ■' 
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However,  the  functions  K (r,t),  are  different  for  the  three  different 


The  bound  obtained  by  Yao  and  Thomas  is 


sinfltl 


where  N,  and  N are  the  number  of  samples  to  the  left  and  to  the  right 


of  the  sample  nearest  to  t.  In  our  notation 


sinittl 


If  we  define  R (r)  as  the  ratio  of  the  bound  obtained  here  to  the 


bound  obtained  by  Yao  and  Thomas,  we  find 


N+m+1/2  N-nri-1/2 


■mi 
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The  ratio.  R.  (r)  for  o < r < 1 is  shown  in  Figure  1 


2.68,  while 


Note  also  that  for  values  of  r > 0.73,  the  bound  on  the  truncation 


error  given  here  is  tighter  than  the  Yao-Thomas  bound.  That  is. 


for  samples  taken  nearly  at  the  Shannon  rate,  the  bound  obtained 


here  is  an  improvement  on  the  Yao-Thomas  bound 


The  bound  given  by  Brown  in 


where  N.  and  N,  are  the  same  as  in  the  Yao-Thomas  bound.  Forming 


the  ratio  of  the  bound  obtained  here  to  the  bound  obtained  by 


Brown,  we  find 


Nfmfl/2  N-tn+1  / 2 
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Ratio  Between  Derived  Truncation  Error  Bound  and 
Yao -Thomas  Bound 


r 


II 1BWIII— II  ■ I W 


-18- 
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Hence,  Che  bound  obtained  here  is  an  improvement  on  the  bound 
obtained  by  Brown  for  all  values  of  the  guard-band  coefficient,  r. 


197  ** 
r-1018 


i 
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g(y)  - 


Jbji 


•in  ^ for  _*  < y < 0 <|b|<  We  want  to 


•in  y 

show  that  g(y)  is  monotone  increasing  on  (-n,n).  Note  that  g(y) 
is  an  odd  function  of  y,  and  hence  it  is  sufficient  to  consider 
o < y < n.  Also,  g(y)  > o for  d < y < « and  g(o)  - 0.  Thus, 
we  only  need  to  consider  o < y < n. 


g (y) 


2 |b I sin  [b|y  cos |b|  y sin  y - cos  y sin  |b|y 


. 2 

sin  y 

- s,in[b  ljr„.  1 2 |b  | cos  |b  | y sin  y - cos  y sin  I b I y 
sin  y [_ 


Case  1:  o < y < j . Then  o < | b |y  < |b  I \ < and  hence. 


sin  |b  | y > o 


cos  | b | y > o 


sin  y > o 


cos  y > o. 


g’(y)  > 


2 lb  I sin  y cos  Ibly  - lb  I y cos  y 


sin  y [_ 
since  sin  Jb  | y < lb  | y. 


Thus, 


g (y)  > °*n  |b  | y 2 8in-^  cos  lb  I y - cos  y | 

sin  y L y J 


_ sin  b ly  L I n a in  2 I.  I 

> — U2U1|b|y  2 “T“  cos  M 

•in  y L 2 


y - cos  y 


> — — l-"-l^  |b  | y cos  jb  | y - cos  y ^ 


> o, 


A 


■[•-  i] 


since  o<|bJy<y<-|  and  cos  y is  monotone  decreasing 
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Case  2:  - < y < «.  Then, 


< -^b|n  ■ < j b|y  <|b|n  < and  therefore 


sin  y > o 


sin|b|y  > o 


' . . sin  b 

g (y)  - — 0 


cos  y < o 
cos|b|y  > o. 


2 ( b I sin  y coslbly  - sin|b|y  cos  y 

sin  y 


Hence,  for  o < y < n,  g(y)>o  and  consequently  g(y)  is  monotone 
increasing  on  -it  < y < «. 

Now  consider 


h(y)  - 


sin  by  cos  b\ 


sin  y 


sin  2by 
2 sin  y 


Note  that  h(y)  is  an  even  function  of  y,  and  therefore  if  h(y)  is 
monotone  Increasing  for  o < y < «,  it  will  be  monotone  decreasing  on 
-»  < y < o.  Also,  if  h(y)  is  monotone  increasing  for  b > o,  it  will 
be  monotone  decreasing  for  b < o.  Thus,  it  is  sufficient  to  consider 
o < y < «,  b > o. 

h '(y)  m 2b  cos  2by  sin  y - sin  2by  cos  y 
2 sin^  y 


Case  Is  o < y < ^ • o < 2by  < ^ . 


4 
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Note  that  h (o)  ■ o,  and  hence  it  is  sufficient  to  consider  o < y < -z 


h (y) 


2b  tan  y - tan  2 fry 

2 tan  v -slHX.  . 

Z Can  y cos  2by 


2 tan  y sln— 7^  , > o for  o < y < ^ , o < 2 by  < |.  Thus,  we  only 

J cos  2 by  2 7 2 


need  consider  2b  tan  y - tan  2 by  , and 


2 b tan  y - tan  2 by  = 2 by 


[tan  y tan  2 by|  . 

y " 2 by  J - °' 


since  is  monotone  increasing  on  o < y < | . Hence, 

h (y)  > o for  o < y < . 


Case  2;  -g  < y < it.  If  o < 2 by  < -g,  then 


cos  2 by  > o sin  y > o 


cos  y < o 


sin  2 by  > o , and  thus 


h' (y)  2 b cos  2 by  sin  y - sin  2 by  cos 

2 sin2y 


> o . 


For  75  < 2 by  < n , -g  < y < n > we  have  cos  2 by  > cos  y , since 


2 by  < y , and  hence. 


h (y)  > C°°  2 b sin  y - sin  2 by 

” 2 sin  y 


2by  cos  y sin 
2 sin^  y y 


sin  2 bi 
2 by 
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sin  2 by  ^ sin  y , ^ , , sin  y , 

But  — g by  — — y — * since  2 by  < y and  — is  monotone  decreasing 

on  o < y < it.  Also,  cos  y < o , and  thus, 

h (y)  > ° > for  \ < y < «. 


We  have  shown  that  h (y)  > o for  o < y < jr  , and  consequently  h(y) 
is  monotone  increasing  on  o < y < n for  b > o . 
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t. 


El 


i &. 


! I 


f ! 1 


Appendix  II 


We  wish  to  evaluate  I I£  - I^l , where  Ig  and  I^  are  given  by  (24)  and 


(25). 

b 


/ sin(2lfrl)y  cos  2 my  dy  = \ J [sin(2N+2rot-l)y  + sin(2N-2m-l)y  dyj 

O n — 


1 J cos(2N+2gn-l)  b - cos(2N+2mfl)  a 

2 1 2N+2hh-1 


cos(2N-2m+l)  b - cqs(2K-2eh-1)  a 

2N-2iih-1 


Also, 

b 


j sin(2N+l)y  sin  2 my  dy  = 


1 I sin(2N-2mfl)b  - sin(2N-2rofl)a 

2 I 2N-2iih1 


sin(2N+2rot-l)b  - sin(2N+2ra-l)a 
2N+2mt-l 


Using  these  integrals  and  simplifying  by  means  of  elementary  trigonometric 
identities,  we  find  that 

slnOh-mi). 


1 ) sin  b it  sin  b 


05 


2 ~ 0 05 

2 cos  75 


2N+2mfl 


+ l^l-cos  bn  cos  booj^cos  (2N+2m+l)B^  + i sin  (2Nt2nrt-l)  i 

2 N+2nri-l 

sin  bn  sin  bo5  sin  (N-nrt-2)  n [je~*  m+j)  tn j 

2 N-2mtl 

[l-  cos  bn  cos  bo^  [cos  (2N-2mfl)p^  + 1 sin  (2N-2mfl)p^J 

2N-2mfl  1 

J 


p I Jl.  J..JJJJ  U.WmiJIJ- 

and 


Forming 


3 eg. 

2 cos  2 
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sin  bn  cos  bca  sin  (tfhnf  l/2)nLle  enj 

2 N+2m+l 


sin  bn  cos  ben  s 


In  (N-mfl/2)n  Cie1^’"*1^ 


2 N- 2m+l 


sin  bn  cos  ben  + cos  bn  sin  b cn 
2 


cos  (2N+2mf  1) P~  + i sin(2N+2m+l)p  cos  (2N-2m+l)p  - 1 sin  (2N-2mfl)p,.  j 

‘ j 


2N+2m+l 


2N-2m+l 


sin  bn  cos  ben  - cos  bn  sin  ben 


cos  (2N+2mfl)P4  + 1 sin  (2N+2m+l)p6  cos  (2N-2m+l)P4  - 1 sin  (2N-2mfl)p6 
2N+2m+l  + 2N-2m+i  t 


I2  " il^*  collecting  terms  and  taking  magnitude,  we  get 


1X3  < 


en_  ) 2N+2nri-l 
2 cos  2 


sin  bnl+2  (1-  cos  bn  cos  ben) 


+ E 

11  o 


sin  bn  cos  ben  + cos  bn  sin  ben 


'IT 


+ ~2 


sin  bn  cos  ben  - cos  bn  sin  ben 


a 1 


2N-2m+l 


jsin  bn  sin  b cn  I +JT  (l-cos  bn  cos  beo) 


IT 


+ |sin  bn  cos  bo>|  (sin  bit  cos  bio  .+  cos  bn  sin  ba> 


♦f  i 


sin  bn  cos  ben  - cos  bn  sin  i 


k“0i 


But 


I 

I 
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2 [ 2 

sin  bn  I = Isin  bn  I Isin  ben  I + cos  lb ao 


If 

r|[s 


1] 


< sin  bit  I sin  ba>  I + cos  bco 


= Isin  bn  sin  bad  + sin  bn  cos  ben 


•Ij 


Isin  bn  cos  ben  + cos  bn  sin  ben  I = sin  b(n-fcn  = sin  lb  I (n+en  ) , 


■M' 


Also, 


and  |sin  bn  cos  ben  - cos  bn  sin  ben  J = | sin  b(n-cn)j  = sin[b|(n-en) 

Or, 

|sin  bn  cos  ben  + cos  bn  sin  ben  | + jsin  bn  cos  ben-  cos  bn  sin  ben 

2 Isin  bn  | cos  b en 


Also, 


2 2 

1 - cos  bn  cos  b en  < 1 - cos  bn=  sin  bn. 


Making  use  of  these  inequalitities , we  find 


I2-U3 


< (1-4-  2)  j sin  bn 


o 0> 

2 cos  j 


1 + 1 


2N-2mfl  2n-2m+l 


f 


I + 


sin  b en  +|  cos  b eo 


which  is  our  desired  result. 
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